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Abstract 

We have carried out 12 CO(J =2-1) and 12 CO(J =3-2) observations at spatial 
resolutions of 1.0-3.8 pc toward the entirety of loops 1 and 2 and part of loop 3 in 
the Galactic center with NANTEN2 and ASTE. These new results revealed detailed 
distributions of the molecular gas and the line intensity ratio of the two transitions, 
i?3_2/2-i- In the three loops, the line intensity ratio R3-2/2-1 is in a range from 
0.1 to 2.5 with a peak at ~ 0.7. The ratio in the disk molecular gas is in a range 
from 0.1 to 1.2 with a peak at 0.4, significantly smaller than that in the Galactic 



center. This supports that the loops are more highly excited than the disk molecular 
gas. An LVG analysis of three transitions, 12 CO J =3-2 and 2-1 and 13 CO J =2-1, 
toward six positions in loops 1 and 2 shows density and temperature are in a range 
10 22 - 10 4 ' 7 cm -3 and 15-100 K or higher, respectively. Three regions extended by 
50-100 pc in the loops tend to have higher excitation conditions as characterized by 
-2/2-1 greater than 1.2. The highest ratio of 2.5 is found in the most developed 
foot points between loops 1 and 2. This is interpreted that the foot points indicate 
strongly shocked conditions as inferred from their large linewidths of 50-100 km s" 1 , 
confirming the suggestion by Torii et al. (2010b). The other two regions outside the 
foot points suggest that the molecular gas is heated up by some additional heating 
mechanisms possibly including magnetic reconnection. A detailed analysis of four 
foot points have shown a U shape, an L shape or a mirrored-L shape in the b-v 
distribution. It is shown that a simple kinematical model which incorporates global 
rotation and expansion of the loops is able to explain these characteristic shapes. 
Key words: radio lines: ISM - ISM: clouds - magnetic loops 

1. Introduction 

The central region of the Galaxy has an outstanding concentration of molecular gas 
within ~300 pc in radius called the Central Molecular Zone (CMZ). The molecular gas shows 
high temperatures and large velocity dispersions, considerably enhanced from those in the 
Galactic disk (e.g., Morris and Serabyn 1996, Glisten & Philipp 2004). In addition, there are 
a few other molecular features having large velocity dispersions of 15 km s -1 or more outside 
the CMZ within the central lkpc, including Clumps 1 and 2 (Bania 1977; Bania et al. 1986 
and Stark et al. 1986) and the I = 5.5° Complex (Bitran et al. 1997). The origin of these high 
temperatures and large velocity dispersions has been a puzzle over a few decades. 

The most outstanding feature, the CMZ, harbors giant molecular clouds including the 
Sgr A and Sgr B2 molecular clouds (e.g., Fukui et al. 1977, Scoville, Solomon and Penzias 1975) 
and, in addition, the CMZ appears to be associated with two extreme velocity features called 
an expanding molecular ring (EMR). These features, the CMZ and EMR, form a parallelogram 
in a longitude-velocity diagram. Theoretical studies by Binney et al. (1991) showed that the 
kinematics of the CMZ and EMR can be reproduced by the non-circular motion driven by a 
stellar bar potential. The model was applied to these features including the CMZ (e.g., Liszt 
2006), whereas the origin of the whole anomalous kinematic features in the central kpc disk 
still remains unsettled (see for a review Morris and Serabyn 1996). 

It has been shown that magnetic field may also play an important role in the Galactic 
center as demonstrated by the Radio Arc (Yusef-Zadeh, Morris and Chance 1984) and non 
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thermal filaments (Yusef-Zadeh, Hewitt and Cotton 2004). These features suggest that the 
magnetic field can be as strong as mG, whereas a large volume with lower density may have 
weaker fields as suggested by low frequency radio observations (LaRosa et al. 2005). A recent 
reanalysis of radio nonthermal radiation indicates a strong field of ~50 //G over ~400 pc scale 
(Crocker et al. 2010). The magnetic pressure may be important in gas dynamics since the 
energy density of the magnetic field is comparable to the kinetic energy density of the violent 
motion. Polarization measurements at sub-mm wavelength indicate that the field direction is 
nearly parallel to the plane (Novak et al. 2003), while vertical non-thermal features indicate 
vertical field directions. The global field configuration yet remains to be clarified. 

Fukui et al. (2006, F06) analyzed the NANTEN survey of the CO (J =1-0) transition 
toward the Galactic center (Mizuno & Fukui 2004) and discovered two huge molecular loops 
in I = 355°-358°. These loops show enhanced velocity dispersions of 50-100 km s -1 and F06 
proposed that they are magnetically floated loops up to heights of 200-300 pc from the plane 
at 8.5 kpc created by the Parker instability (Parker 1966). F06 presented a picture that the 
Parker instability formed two loops and the gas falling down along the field lines accumulated as 
two "foot points" of molecular gas toward each end of the loops, with an estimate of magnetic 
field strength of 150 /xG and an Alfvenic speed around 24 km s -1 . Subsequently, Fujishita 
et al. (2009) discovered another molecular loop, loop 3, and argued that it also represents a 
magnetic flotation loops. Torii et al. (2010a and 2010b, hereafter TlOa and TlOb) made follow 
up observational studies of loops 1 and 2 and a detailed analysis of muti- J CO transitions in 
TlOb derived temperature of 30-100 K or higher and density of 10 3 cm -3 in the foot points and 
TlOb argued for shock heating in the foot points. Numerical simulations of magnetic loops were 
also made in a global kpc scale disk (Machida et al. 2009) and in a 2-D local loop (Takahashi 
et al. 2009), offering a theoretical basis of the magnetic loops. 

These recent studies of Galactic magnetic flotation have promoted our understanding of 
the properties of the loops in the Galactic center since 2006. It is important at this stage to 
reveal the temperature and density distributions not only toward foot points but also over the 
entire loops. For this purpose we have carried out new high resolution molecular observations 
by using the NANTEN2 4m sub-mm telescope and the ASTE 10 m sub-mm telescope in the 
12 CO(J=2-l), 13 CO(J=2-l), and 12 CO(J=3-2) transitions for loops 1 and 2 and part of loop 
3. Distance is taken as 8.5 kpc (IAU 1985 standard; Kerr and Lynden-Bell 1986). This paper 
is organized as follows; the observations are presented in section 2 and the results in section 3. 
Details of the data analysis are shown in section 4, and discussion is given in section 5. We 
conclude the paper in section 6. 
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2. Observations 

2.1. NANTEN 2 

Observations of the 12 CO(J =2-1) and 13 CO(J =2-1) transitions were made with the 
NANTEN2 4m sub-mm telescope of Nagoya University at Pampa la Bola, Chile, at an altitude 
of 4865 m, in February, October-November 2008 and June 2009. The half-power beam width 
(HPBW) was ~90 arcsec at 230.538 GHz the frequency of 12 CO(J =2-1) and at 220.398681 
GHz the frequency of 13 CO(J =2-1) as measured by observing the Jupiter. We observed the 
entirety of loops 1 and 2, and two parts of Loop 3 (the west foot point and part of the loop 
top) in 12 CO(J =2-1) by using the OTF (= on the fly) mapping technique and six positions 
in loops 1 and 2 in 13 CO( J =2-1) using a position switching mode. NANTEN2 was equipped 
with a 230 GHz SIS receiver that provided a typical system temperature of 150-230 K in 220- 
230 GHz in the double-side band, including the atmosphere toward the zenith. The intensity 
calibration was made using the chopper-wheel method. For checking the whole system and 
absolute intensity calibration, M17SW [R.A.(1950) = 18 h 17 m 30?0, Dec.(1950) = -16°13'6."0] 
or p-Ophiuchus [R.A.(1950) = 16 h 29 m 20?8, Dec.(1950) = -24°22'13"5] were observed every 
2-3 hours. 

We observed Ori-KL [R. A. (1950.0) = 8 h 45 m 14f8, Dec. (1950.0) = 13°30'40"0 or I 
= 213.°6695, b = 31.°7922] and the Rosetta molecular cloud [R. A. (2000.0) = 6 h 34 m 6f9, 
Dec.(2000.0) = 4°25'37. // 2 or / = 207.°0154, b = -l.°8228] with NANTEN2 for the absolute 
intensity calibration. We compared the spectra obtained at Ori-KL and Rosette with those ob- 
tained with KOSMA (Schneider et al. 1998) and the 60 cm telescope of the University of Tokyo 
(Nakajima et al. 2007), respectively, where all the spectra were smoothed to a beam size of 
130" in FWHM. We then obtained the main beam efficiency of 85 % of NANTEN2 at 230 GHz. 
This efficiency was used to convert the observed value to the main beam temperature. The rela- 
tive uncertainties in the calibration are estimated to be better than —5% both in 12 CO(J=2-l) 
and 13 CO( J =2-1) transitions as estimated from the daily variation in the observed brightness 
temperature of M17SW [R.A.(1950) = 18 h 17 m 30f0, Dec.(1950) = -16°13'6"0 or I = 15.°0293, 
b = -0.°6722] or /^Ophiuchus [R. A. (1950) = 16 h 29 m 20?8, Dec. (1950) = -24°22'13."5 or / = 
354.°0554, b = 15.°9400]. 

The spectrometer was an acousto-optical spectrometer with 2048 channels. The fre- 
quency coverage was 250 MHz, corresponding to a velocity coverage of 390 km s™ 1 and a 
velocity resolution of 0.18 km s _1 at 230 GHz. The telescope pointing was measured to be 
accurate to within 10" by radio observations of the Jupiter. The final rms noise level in T m b is 
~0.23 K in 12 CO(J =2-1) and -0.04 K in 13 CO(J =2-1) at 0.18 km s^ 1 resolution. 
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2.2. ASTE 

Observations of the 12 CO (J =3-2) line were made using the Atacama Submillimeter 
Telescope Experiment (ASTE) 10m sub-mm telescope of NAOJ (National Astronomical 
Observatory of Japan) at Pampa la Bola, Chile, at an altitude of 4865m (Ezawa et al. 
2004, Ezawa et al. 2008) in August-September 2006 and April- June 2008. The HPBW 
was 22" at 345.79599 GHz, the frequency of 12 CO(J =3-2). We observed the entirety of 
loops 1 and 2 by using a position switching mode in 2006 and an OTF mapping technique 
in 2008. ASTE was equipped with a 345 GHz SIS receiver, SC 345, that provided a typ- 
ical system temperature of 190-300 K in 2006 and CATS 345 that provided a typical sys- 
tem temperature of ~150 K in 2008, in the double-side band including the atmosphere to- 
ward the zenith. The beam efficiency was measured to be 60 % for 345 GHz (ASTE web 
site; http:/ /www. nro.nao.ac.jp/~aste/prop08/aste_status. html). The intensity calibration was 
made using the chopper-wheel method. M17SW [R.A.(1950.0) = 18 h 17 m 30?0, Dec.(1950.0) 
= -16°13'6"0 or / = 15.°0293, b = -0.°6722] or NGC6334 [R.A.(1950.0) = 17 h 17 m 32f2, 
Dec.(1950.0) = -35°44'0"4 or I = 351M172, b = -0.°6465] was observed every 2-3 hours 
and the absolute intensity scale was established by comparing the intensity of M17SW between 
ASTE and CSO (Wang et al. 1994). The spectrometer consisted of four digital "XF-type" 
spectrometers with 2048 channels with a frequency coverage of 512 MHz, corresponding to a 
velocity coverage of 450 km s _1 and a velocity resolution of 0.43 km s _1 at 345 GHz. The 
telescope pointing was measured to be accurate to within 2" by radio observations of the 
Jupiter, and W-Aql [R.A.(2000.0) = 19 h 15 m 23?l, Dec.(2000.0) =-7°2'49."8 or I = 29.°7224, b = 
-9.°0720] in 2006, and RAFGL5379 [R.A.(2000.0) = 17 h 44 m 23?9, Dec.(2000.0) = -31°55'39"1 
or / = 357.°6555, b = -l.°9358] and IRAS 16594-4656 [R.A.(2000.0) =17 h 03 m 10?0, Dec.(2000.0) 
= -47°00'27"6 or I = 340.°7232, b = -3.°8346] in 2008. The final rms noise level was ~0.33 K 
in T mb in 12 CO( J =3-2) at 0.43 km s" 1 resolution. 

3. Results 

3.1. Distributions of the molecular gas 

Figure 1 shows the integrated intensity distributions of 12 CO (J =2-1). The velocity 
range is —180 — 90 km s _1 including the main body of loop 1 in Figure 1(a) and —90 — 40 km 
s" 1 including the main body of loop 2 in Figure 1(b) (F06). Figures 2(a) and (b) show the 
intensity distributions in the same two velocity ranges in 12 CO(J =3-2). Figures 3(a) and (b) 
show loop 3 for a velocity range 30-160 km s _1 in the integrated intensities of 12 CO(J =2-1) 
and 12 CO(J=3-2), respectively. The solid thick contours in Figures 1-3 show the lowest 3cr 
NANTEN 12 CO(J=1-0) contour, 7.1 K km s" 1 (TlOa). 

The 12 CO(J =2-1) and 12 CO(J =3-2) transitions generally show similar distributions 
to 12 CO(J=1-0), while the 12 CO(J=2-l) and 12 CO(J=3-2) intensities seem weaker than the 
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12 CO(J =1-0) intensity especially in part of the tops of loops 1 and 2. The 12 CO(J =3-2) 
distribution allows us to find detailed distributions at 1 pc resolution. In each panel of Figures 
1 and 3, the three regions of the foot points of loops 1 and 2 and the western foot point of loop 
3 are clearly seen; the eastern foot point of loop 1 is located in (7,fc)~(357.°2-357.°6, 0.°4-0.°9) 
(Figures 1(a) and 2(a)), the western foot point of loop 1 and the eastern foot point of loop 2 
in (Z,&)~(355.°8-356.°3, 0.°6-l.°3) (Figures 1(b) and 2(b)), the western foot point of loop 2 in 
(/,6)~(355.°2-355.°6, 0. 8-l.°2) (Figures 1(b) and 2(b)), and the western foot point of loop 3 in 
(1, 6)~(355.°3-355.°8, 0.°5-0.°8) (Figures 3(a) and (b)). The present observations do not cover 
the eastern foot point of loop 3. We present detailed velocity channel distributions every 10 km 
s _1 of the two transitions in Figures A and B in the Appendix 1 of Supporting Online Material. 
We shall denote the observed area by eight regions for discussion as shown in Figure 4; Regions 
A, B and C include loop 1, Regions C, D and E loop 2, and Regions F and G loop 3. 

3.2. Line intensity ratios 

The line intensity ratio 12 CO(J =3-2)/ 12 CO(J =2-1) (hereafter # 3 _ 2 /2-i) is shown for 
loops 1 and 2 and for loop 3 in Figures 5(a) and (b), and Figure 6, respectively, for the same 
velocity range as in Figures 1 and 2. The 12 CO( J =3-2) data is smoothed to the same Gaussian 
beam of 90" FWHM of the 12 CO(J =2-1) data. The ratio is calculated for positions where 
the two line intensities are larger than the 3a level. Relatively high values of #3-2/2-1 around 
0.7 to 1.3 are seen in most of the loop 1 foot points in addition to the two regions toward the 
western inner side of the loop 1 top around (Z,6)~(356°.5, l.°2) and the western loop 2 top 
around (/,6)~(355.°3, l.°5). We further note that two small regions toward (Z,6)~(356.°0, l.°5) 
and (355.°9, l.°9) show high values of 1.3 in #3-2/2-1- The eastern foot point of loop 1 is not so 
high in #3-2/2-1- Detailed distributions of #3-2/2-1 every 10 km s _1 are presented in Figure C 
in Appendix 1 of Supporting Online Material 1 . In Figure 5 we also indicate six areas of high 
line intensity ratio, P, Q, R, S, T and W as listed in Table 1. 

In Figure 7 we show intensity-weighted frequency distributions of #3-2/2-1 for the ve- 
locity range of the loops, —140 — 40 km s -1 , in white and the three intervals of local features 
in —40 — 20 km s _1 in orange, —20-0 km s _1 in red, and 0-20 km s" 1 in green for the data 
enclosed by the broken lines in Figure 5. The grey area in Figure 7 indicates #3-2/2-1 i n the 
velocity interval of the loops, where #3-2/2-1 ranges from 0.3 to 1.2 at a 10 % level of the 
maximum value at 0.7. In the three local velocity intervals from —40-20 km s _1 , #3-2/2-1 
histograms have peaks at 0.4-0.5. This difference in # 3 _ 2 /2-i suggests that the gas in the loops 
have more highly excited conditions than in the foreground local features, whereas the lower 
velocity gas may include some of the loop components. 

Figure 8 shows intensity weighted frequency distributions of #3-2/2-1 i n the seven in- 

1 The fits files of loops 1 and 2 in 12 CO(J=2-l) and 12 CO(J=3-2) can be taken from 
http : / / www. a. phys . nagoya-u. ac . j p/ ~ torii/loop/cg 
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dividual regions in Figure 4. Figure 8(a) is for loops 1 and 2 and Figure 8(b) for loop 3. In 
Figure 8(a) Region C with the most intense foot points show a wing-like distribution reaching 
2.5 in i?3_2/2-i, confirming the high excitation states found by TlOb. Region B shows a similar 
distribution up to 2.2 and Region E up to 1.7. Regions B, C and E have peaks at 0.5-0.7 
while Region D is unique in that its peak is high, around 1.3. Region A shows less values of 
#3-2/2-1, whereas the observed region in 12 CO(J=2-l) in Region A does not cover the southern 
part of loop 1 toward (Z,&)~(357.°5-358.°0, 0.°0-0.°4). The high values of #3-2/2-1 suggest high 
molecular gas temperature and/or density and indicate that four Regions, B, C, D and E are 
in higher-excitation states. The intensity weighted frequency distributions of #3-2/2-1 m 1°°P 
3 are shown in Figure 8(b). Generally, #3-2/2-1 is less than 1.2 peaked at around 0.6 and the 
molecular gas in loop 3 is of lower excitation than those in loops 1 and 2. 

These results indicate that Regions C and D exhibit highest excitation conditions in 
the observed area and Region B second highest. Region C was already analyzed by TlOb. It 
is notable that the other foot points, Regions A and G, show lower excitation conditions and 
Region E has only slightly higher excitation conditions than Regions A and G (see the #3-2/2-1 
distribution in Figure 8). So, it is not true that all the foot points show very high excitation 
conditions. None the less, we note that their #3-2/2-1 distributions peaked at around 0.6 are 
still of higher excitation than the disk clouds that have #3-2/2-1 peaked at around 0.4-0.5 
(Figure 7), and that the molecular gas in the loops show higher excitations than in the disk. 

3.3. Details of five regions 

We describe more detailed properties of the five regions, A, C, D, E and G, in the 
following. 

Region A (Loop 1 eastern foot point) 

Figure 9(a) shows velocity channel distributions of the 12 CO(J =3-2) intensity and 
#3-2/2-1- We find filamentary distributions along the loop in —140 — 130 km s _1 in the 
12 CO(J =3-2) while the 12 CO(J =2-1) observations cover only the northern part of the fila- 
ment. The area around (/,6)~(357.°5, 0.°8) shows that a CO(J=1-0) clump in TlOa is resolved 
into a V-shaped feature and has no point whose #3-2/2-1 is larger than 1.2. Figure 9(b) shows 
velocity-position distributions along the filament. We define a Y axis along the filament in 
an area of (/,6)~(357.°5-357.°9, 0.°0-0.°4) by taking a least-squares fitting with a straight line 
weighted by the integrated intensity. An X axis is taken vertical to the Y axis as shown in 
Figure 9(b) left. Figure 9(b) right shows the integrated intensity every 0.1° along the X axis 
in velocity-position distributions. Below Y = 0.4°, we see only narrow features but above Y = 
0.4° we find the features nos.l and 2 has broad velocity distributions from —180 to —130 km 
s _1 as indicated by dotted lines. The overall shape is like a "mirrored L" in —160 — 120 km s _1 
and Y = 0.°8-l.°l as indicated by dash-dotted line. We discuss these features in Section 5.2. 
Region C (Loop 1 western foot point) 
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Figure 10(a) shows velocity channel distributions of the 12 CO(J =3-2) intensity and 
#3-2/2-1- We find two regions of high 12 CO( J =3-2) intensity; one of them in a velocity range 

—100 80 km s" 1 in I = 356.2° and b = 0.85° shows a strong intensity gradient toward the 

plane and the other in —60 - 50 km s _1 in / = 356.2° and b = 0.7° shows a protrusion toward 
the plane. The gas in a velocity interval —60-10 km s _1 in (Z,6)~(356.°2, 0.°9) shows #3-2/2-1 
greater than 1.7 and this corresponds to the point inside the U shape discovered by TlOb. This 
feature with high #3-2/2-1 lead TlOb to suggest that the foot point of loop 1 is extended to 
— 10 km s" 1 instead of —50 km s" 1 as originally suggested (F06) in a l-v distribution (Figure 
2 in F06). There is a hot spot named Q in an interval —100 — 80 km s _1 in (/,6)~(356.°1, l.°5) 
which shows a high #3-2/2-1 ratio around 1.7. This feature has been discovered in the present 
work, since TlOb does not cover b greater than l.°2. 

Figure 10(b) shows b-v distributions of the foot point. TlOb discovered a U shape in the 
integration range from I ~ 356.°15 to 356.°21 as indicated by dash-dotted lines and found four 
broad emission features, nos. 3, 4, 5 and 6, as shown by dotted lines in Figure 10(b). Among 
them the feature at I = 356.° 15-356.° 17 (no. 5) shows temperature higher than 100 K (see section 
4.2 of TlOb). In a range from / = 356.°15 to 356.°19 in Figure 10b a protrusion is indicated by 
bold arrows, showing a velocity dispersion of 30 km s -1 . 
Region D (Loop 2 top) 

Figure 11 shows Region D in 12 CO( J =3-2) in a velocity range —90 — 60 km s _1 . We find 
the helical distribution noted by TlOa toward (Z,6)~(355.°72, l.°80) and filamentary features, 
some of which are elongated normal to the loop direction toward (/,6)~(355.°88, l.°92), (355.°54, 
l.°80) and (355.°40, l.°76). We find these features show higher excitation conditions toward a 
hot spot S in the helical distribution and toward the filamentary features in T and W in the 
lower panel of Figure 5. 
Region E (Loop 2 western foot point) 

Figure 12(a) shows velocity channel distributions of Region E in 12 CO(J =3-2) and 
#3-2/2-1; where we see a protrusion in —90 — 70 km s" 1 toward (/,6)~(355.°4, 0.°7). Figure 12(a) 
indicates that Region E shows #3-2/2-1 ratio less than ~1.7 and #3-2/2-1 of this protrusion is 
highest in Region E. Figure 12(b) shows a b-v distribution of the region. We see an L-shaped 
feature in I ~355.°42-355.°48 as indicated by a dash-dotted line. We also see three regions with 
broad linewidths, nos 7-9, indicated by dotted lines and the protrusion noted by bold arrows 
in Figure 12(a). This distribution has a velocity dispersion of ~10 km s _1 . 
Region G (Loop 3 western foot point) 

Figure 13(a) shows a velocity channel distributions of Region G in 12 CO (J =3-2) and 
#3-2/2-1- The gas comes closer to the left with the increase of velocity from 50 to 100 km s _1 , 
i.e., showing a velocity gradient of 2.0 km s _1 pc _1 . We see two concentrations of gas toward 
v -60-70 km s" 1 and (Z,6)~(355.°5, 0.°7) and v -80-90 km s" 1 and (Z,6)~(355.°6, 0°6) where 
#3-2/2-1 is enhanced. Like Region A we find no gas with #3-2/2-1 larger than 1.7. Figure 13(b) 
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shows a velocity-position diagram in the integration range over I ~355.°46-355.°53. We see a 
mirrored-L shape indicated by a dash-dotted line. We find five regions with broad linewidths, 
nos. 10-14, as noted by dotted lines but we do not see small protrusions extended in the north 
to south like that in Region C or E. 

4. Data analysis 

Figure 14 shows a plot of the line intensity ratio 
R^2/2-icalculatedbytheLVGapproximation{GoldreichandKwan 1974) as a function of 
temperature of 10-200 K for density from 10 2 cm -3 to 10 4 cm -3 . In these LVG calculations 
(Goldreich and Kwan 1974), i? 3 _ 2 /2-i ranges from ~0.1 to ~1.0 and is a monotonic function 
of density and temperature. -R3-2/2-1 becomes greater than 0.7 for density higher than ~ 10 3 
cm -3 and temperature higher than ~70 K. These conditions appear to characterize the high 
excitation gas in the Galactic center. The typical cloud condition in the Galactic disk, 10 2 -10 3 
cm -3 and 10-20 K, corresponds to #3-2/2-1 less than 0.5. These properties are consistent with 
the trend of i? 3 _2/2-i distribution in Figure 7. This confirms that i? 3 _ 2 /2-i is an indicator 
of the excitation condition, a combination of density and temperature, while one cannot 
determine the parameters uniquely with these two transitions only. 

Figure 5 and Table 2 show six positions where pointed observations were made in the 
13 CO( J =2-1) transition with NANTEN2; they are a in Region A, b in Region B, cl and c2 in 
Region C, d in Region D, and e in Region E. Figure 15 shows the line profiles of the six positions, 
including 12 CO(J=2-l), 13 CO(J =2-1) and 12 CO(J=3-2). The 13 CO(J =2-1) emission was 
significantly detected in all the positions. Table 3 lists the ratios in the integrated intensity 
every 10 km s^ 1 of 12 CO( J =3-2) to 12 CO( J =2-1) (#3-2/2-1) and those of 13 CO( J =2-1) to 
12 CO(J =2-1) (i?i3/i2)- In order to estimate the physical parameters of these positions, we 
used an LVG analysis (Goldreich and Kwan 1974) for the 12 CO(J=2-l), 12 CO(J =3-2) and 
13 CO(J=2-l) transitions. 

4-1. LVG analysis 

In order to understand the physical conditions in molecular clouds we need to solve 
simultaneously the equations of line radiation transfer and statistical equilibrium of the molec- 
ular energy levels by taking into account radiative and collisional processes. In a large velocity 
gradient (LVG) analysis we assume a sphere with a uniform velocity gradient and uniform 
temperature and density. A photon emitted at any position can escape from the cloud at an 
escape probability (3 = (1 — exp(— r))/r, a function of line optical depth r, that is multiplied to 
the Einstein coefficients to correct for photon escaping due to the optical depth effect (Castor 
1970). We then estimate the line intensity of a particular transition for kinetic temperature, 
density and molecular abundance X(CO) / (dv / dr) , where X(CO) is the ratio of CO to H 2 and 
(dv/dr) the velocity gradient, by combining statistical equilibrium equations of CO rotational 
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levels, with Einstein coefficients and collisional excitation rates (Schoier et al. 2005). Here we 
assume that X{CO) / {dv / dr) is 1.0 x 10~ 5 (TlOb), [ 12 C]/[ 13 C]~ 53 (Wilson & Rood 1994; see 
also Riquelme et al. 2010) and [ 12 CO]/[H 2 ] ~ 10~ 4 (e.g., Lis et al. 2001). 

4-2. Results 

We have calculated the line intensities over a range of physical parameters for the line 
intensity values listed in Table 3 and calculated x 2 distributions between calculated ratios and 
observed ratios. The degree of freedom is 1 for three transitions. Figure 16 shows the x 2 
distributions of 95% confidence level for the six positions (x 2 = 3.84) where a cross indicates 
the point for minimum x 2 ■ Table 4 lists the range of physical parameters, i.e., maximum and 
minimum values that can be realized at 5% probability. Figure 17 shows the results with line 
profiles in Table 4. Typically speaking, density is ~ 10 3 cm -3 and temperature is in a range 
of 10-50 K. The positions cl and d show highest temperatures of ~100 K or even more. This 
is generally consistent with the distribution of i? 3 _ 2 /2-i in Figure 8(a) which indicates highest 
ratios in Regions C and D. 

The present six positions were not selected based on i?3_2/2-i an d have lower R3-2/2-1 
of 0.3-1.5, representing the whole loops, whereas Figure 8 indicates highest i?3_2/2-i reaches 
2.5. We note that the positions with i? 3 _ 2 /2-i larger than 1.0 show temperatures even higher 
than ~ 100 K at the highest density 10 3 ' 5 -10 4 cm 3 in the LVG results (Table 4). Considering 
these and the previous results in the loop 1 foot point by TlOa, we summarize that the loops 
generally have temperatures around 10-50 K with highest spots of above 100 K and the typical 
density is 10 3 cm -3 . These are significantly different high excitation conditions from those of 
molecular gas in the disk that has typical temperature of 10 K when no local heat source like 
high-mass stars is present (e.g., see Mizuno et al. 1995 for the Taurus low-mass star forming 
region) . 

5. Discussions 

5.1. Previous studies 

The present 12 CO (J =3-2) and 12 CO (J =2-1) study of the loops in the Galactic center 
has shown that the loops have higher excitation conditions of molecular transitions which are 
not typical of the gas in the disk. In this section we shall discuss the individual highly-excited 
regions in the loops. 

TlOa carried out a detailed observational study of loops 1 and 2 in the CO J =1-0 
transition at 10 pc resolution and discovered helical distributions similar to the solar helical 
loops, a sign of magnetic instability. TlOa found counterparts to loops 1 and 2 in the Hi gas 
and far infrared dust features and derived new estimates of the total mass and kinetic energy 
of the loops, ~ 2 x 10 6 M Q and ~ 10 52 ergs, respectively. Fujishita et al. (2009) discovered 
a new loop named loop 3 in the positive- velocity range in / ~ 354°-359°. F06 and Fujishita 
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et al. (2009) argued that loops are created by the Parker instability. More recently, TlOb 
analyzed multiple- J CO transitions including CO (J =1-0, J =3-2 and J =4-3) in the foot 
points between loops 1 and 2 and derived temperatures of 30-100 K or higher, suggesting that 
shock heating is responsible for raising the temperature. TlOb also discovered that the foot 
points show a U shape in a position-velocity diagram with a few compact and broad features 
inside the U shape. They argued that the U shape reflects the magnetic field morphology as 
predicted by MHD numerical simulations and the broad features may result from magnetic 
reconnection between the two sides of the U shape that have anti-parallel field direction. 

Theoretical studies were made to simulate formation of the magnetically floated loops. 
Matsumoto et al. (1988) made two dimensional magneto-hydro-dynamical numerical simula- 
tions of the process in the Galactic plane and demonstrated formation of the foot points with 
shock fronts. Machida et al. (2009) made three dimensional simulations of a kpc-scale global 
nuclear disk and Takahashi et al. (2009) made two dimensional simulations in the disk and 
showed that top heavy loops can be generated as is consistent with the observed molecular 
loops. They argued that the mechanism is similar to the formation of the solar dark filament 
seen in H-alpha and their Figure 10 demonstrated that the down flow of gas forms shock, 
heating a large volume in the foot points of a loop. 

5.2. Physical properties of the three loops 

We note that there is no stellar heat source detected in these regions and any type of 
local stellar heating is not applicable here. In all these regions there are no IRAS point source 
candidates for protostars (Takeuchi et al. 2010) and no radio continuum sources at la noise level 
(0.1 Jy beam -1 ) in the radio data at 8.35 GHz in the First Galactic Plane data (GPA, Langston 
et al. 2000). This 8.35 GHz flux density corresponds to 10 46 FUV photons s -1 , equivalent to a 
single Bl star of zero age main sequence (Kurtz et al. 1994, Panagia 1973, Mezger et al. 1967). 
The present LVG analysis suggests that the density is likely in a range of (1-3) xlO 3 cm -3 and 
high values of -R3-2/2-1 indicate higher temperature rather than even higher density. Therefore, 
we require some other mechanism to explain the high excitation condition. 

We summarize four possible mechanisms for the high excitation conditions already dis- 
cussed in the preceding works (Matsumoto et al. 1988, Takahashi et al. 2009, TlOb); 

a) The gas in the loops falls down to the plane at Alfven speed of a few 10 km s~ l to cause 

shocks that compress and heat the gas (Matsumoto et al. 1988). This type of heating is 
supported by a multi-line analysis of the foot points between loops 1 and 2 which derived 
significantly high temperature of 30-100 K or more (TlOb). 

b) Large-scale numerical simulations of the loops show that the gas falling down to the foot 

point bounces to form a spur moving upward vertical to the plane up to an altitude of 
several H (Matsumoto et al. 1988). The upward motion is also in the order of the Alfven 
speed and can cause shock heating in the spur. Such a spur is suggested in between loops 
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1 and 2 (TlOa). 

c) The magnetic reconnection where field lines become contacted in opposite directions can 

release thermal energy that may accelerate gas in the surroundings by pressure gradient. 
This is suggested to explain the broad features within the U shaped foot point in a latitude- 
velocity diagram in TlOb. 

d) The rising motion inside a loop by continuously floating gas from the plane is suggested by 

Takahashi et al. (2009). This motion can cause shock fronts inside the loop. 

The high temperature of foot point Q are explained by mechanism a), as shown in TlOb. 
We shall discuss the other features below. If we assume a typical molecular clump of 2.5 x 10 4 
M Q having density 10 3,5 cm" 3 and temperature 50 K, the molecular cooling rate is estimated 
to be 4.3 x 10 35 erg s _1 (Goldsmith and Langer 1978). Magnetic energy stored in the clump is 
7 x 10 50 erg for 150 /iG. For a hot spot like clumps in R and S with a 10 pc scale, magnetic 
reconnection can release energy of 10 37 erg s _1 for a typical time scale of 1 Myr. The kinetic 
energy of the clump is estimated to be ~ 2 x 10 50 erg for motion at the Alfvenic speed, similar 
to that of the reconnection energy. We thus infer that both gravitational and reconnection 
energies may play a role to heat up the molecular clump by considering uncertainties in the 
above estimate including the unknown efficiency of kinetic energy conversion in reconnection. 

Region A shows apparently weak gas concentration and a low velocity dispersion, sug- 
gesting that its shock may be weaker than in the other foot points. Its small angle to the plane, 
~45 degrees, may explain the weak shock. Fujishita et al. (2009) suggested that loop 3 may 
be in an earlier evolutionary stage than loops 1 and 2. This leads us to suggest that Region G 
is not well developed yet as a foot point, explaining its weak shock. For Region E TlOa noted 
that a mini loop-like distribution is seen above the foot point. We may speculate that the mini 
loop may help to reduce the falling velocity, leading to a softer shock with less excitation. To 
summarize, Region C shows the largest velocity dispersion and is consistent with the highest 
excitation conditions among all the foot points. 

In the followings we shall discuss the other individual high excitation areas. 

i) The inside of loop 1 (P) 

P is located in the south of Region B where the inner part of loop 1 is highly excited. We 
suggest that this local heating may be due to shock on loop 1 by mechanism d) rising gas 
inside the loop. On the solar surface there are a number of loops instead of a single loop 
(Pariat et al. 2004). If we assume that loop 1 also consists of numerous unresolved loops 
at present, the subsequent loops rising at the Alfven speed may cause heating in the lower 
part of loop 1 via shocks. This is consistent with that the layer of high R3-2/2-1 is thin. 
We assume that this high excitation region has size of 75 pc x 15 pc x 15 pc, and also 
assume density of 10 3,5 cm -3 and temperature of 50 K since R3-2/2-1 exceeds 0.7. The 
cooling power of the region is 5.8 x 10 37 erg s _1 and the cooling time scale is ~ 2.0 x 10 4 yr. 
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This suggests that heating is very recent. We shall estimate the possible energy deposit 
by the rising H I gas. The H I column density is estimated by 

iV(Hl) = 1.82 x 10 18 x W(Rl) (cm' 2 ) = 1.82 x 10 18 x TAV (cm" 2 ) (1) 

where py(Hl)(K km s _1 ) is the Hi integrated intensity, T(K) is the main beam brightness 
temperature of Hi and AV(km s" 1 ) is the half power line width. The mass of the Hi 
weak feature is estimated to be M(Hl) ~ 1.0 x 10 5 M Q for a size of 75 pc by 75 pc, in —180 
- -160 km s _1 . If the gas is moving at the Alfven speed, 24 km s _1 , the kinetic power is 
estimated to be ~ 10 39 erg s _1 , much higher than the cooling energy in a given time scale. 

ii) Spur between loops (R) 

R is part of the spur in Region C between loops 1 and 2 (TlOb) and mechanism b) is 
applicable here. For a cross section of 100 pc 2 of the clump, the gas flowing upward at 10 
km s _1 has a kinetic power of 10 36 erg s _1 . This is in the order of the cooling power shown 
above. 

iii) Loop 2 tops (S, T and W) 

In Region D, we see three high excitation features in v ~ —90 — 60 km s" 1 ; one of them 
is named S toward the helical distribution in (l,b) ~(355.°88, l.°92) (TlOa) and the others 
named T and W, the filamentary distributions toward (l,b) ~(355.°54, l.°80 ) and (355.°40, 
l.°76) (Figure 11). These high excitation features are high in altitude and mechanism c) 
seems to be most likely. The helical distribution is able to create opposite field directions 
that is required for reconnection and the filamentary features in T and W may also include 
helical distributions yet unresolved. 

5. 3. Other properties 

5.3.1. The U-shape distributions 

We presented the spatial and velocity distributions in 12 CO(J =3-2) for the four foot 
points in Figure 18. It is recognized that the four foot points show characteristic features in the 
b-v diagrams; i.e., a U shape, or a similar L or mirrored-L shape is commonly seen in the foot 
points though their spatial distributions vary from one to another. These shapes are indicated 
by a dash-dotted line in Figures 9, 10, 12 and 13 and are summarized by a schematic drawing in 
Figure 19. Figure 19(a) shows the three loops projected in the nuclear disk (TlOa, Fujishita et 
al. 2009) and Figure 19(b) a b-v diagram of the four foot points. The two foot points nearer to 
the sun, Regions A and G, have concentration of gas on the larger velocity side and part of the 
gas in the bottom of the U shape tends to extend to the lower velocity side. On the other hand, 
the gas in the two foot points, Regions C and E, is concentrated on the lower velocity side. 
Their apparent velocity extents are ~ 20 km s _1 for Region A, ^60-80 km s -1 for Region C, 
~35 km s -1 for Region E, and ~30 km s _1 for Region G. It is notable that Region C shows the 
largest velocity dispersion, whereas the other three have smaller velocity dispersions of 20-35 
km s . 
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The gas in Region C in general shows high values in R3-2/2-1 over 1.0 (Figure 10a). 
The U-shaped distribution was discovered by TlOb and these authors suggested that the U 
shape may be formed in the magnetic flotation by the down flowing gas on both sides of a loop. 
This interpretation is qualitatively applicable to all the foot points. The velocity dispersions of 
20-35 km s _1 of the three foot points except for Region C are consistent with the Alfven speed 
24 km s~ x estimated in F06. Quantitatively, however, Region C shows a much larger velocity 
dispersion and requires some additional velocity separation as noted by TlOb. The foot points 
of loops 1 and 2 for instance may be separated in space, causing a larger velocity separation 
due to different rotation and/or expansion velocity. 

We here attempt to model the basic features in the foot points by adopting the expansion 
and rotation derived by TlOa and the falling motion of 2-D numerical simulations (Takahashi et 
al. 2009). Details are given in Appendix 2 of Supporting Online Material. The model indicates 
that the kinematics is basically dominated by the expansion and rotation of 50-150 km s _1 
and that the foot point shows U or L/mirrored-L shapes as observed in b-v diagrams (Figure 
19). The model therefore qualitatively reproduces the observations, whereas qualitative details 
like velocity spans are not well fit. This is actually expected because the 2-D model is not 
finely tuned to fit the observations. For instance the latitude of the foot points is too low as 
compared with the observations and the velocity widths are not taken into account. We also 
note that assumptions adopted in TlOa and Fujishita et al. (2009) are too simple, allowing 
only a first-order approximation at best. Nevertheless, we regard the basic shapes like U and 
L/mirrored-L support that the idea of the model is in the right direction. 
5.3.2. Compact broad features 

We shall discuss the broad features that appear to connect two velocity features in U 
shapes of the foot points (Figures 9, 10, 12 and 13 shown by dotted lines). We numbered them 
as 1-14 and list their integrated intensity and line intensity ratios in Table 5. TlOb noted one 
of the features no. 5 in Region C for the first time and more than ten similar features indicate 
that the compact broad features are common in the foot points. 

All the regions with high excitation conditions do not have heating sources (See Section 
5.2). The total gas mass with large velocity dispersions in Region C is estimated to be ~ 1.0 x 10 3 
M for no.5 (TlOb, by using the Mopra 12 CO(J =1-0) data with an X factor of ~ 0.7 x 10 20 
(cm -2 / (K km s _1 ), TlOa) and the total kinetic energy is estimated to be 4 x 10 48 erg for velocity 
dispersion ~ 20 km s" 1 . Magnetic reconnection between field lines in the opposite directions 
may have occured in the U shape as discussed by TlOb. Assuming reconnection is working, we 
shall compare the energy of the magnetic field and the kinetic energy of the gas with a large 
velocity dispersions in Region C. The total magnetic energy is estimated to be 4.0 x 10 50 erg 
(TlOb). If the magnetic energy is converted into this kinetic energy, the conversion efficiency 
is 2 %, indicating that the magnetic reconnection is a viable explanation. 
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5. 3. 3. Protrusions 

Figures 10 and 12 show Regions C and E have protrusions as indicated by a bold arrow. 
These have velocity spans of 20-40 km s _1 over 20 pc and are elongated nearly vertical to 
the plane. In the protrusions, the ratio #3-2/2-1 is around ~0.7, higher than that in the 
surroundings of ~0.4. We speculate that the origin of these protrusions is due to downward 
force either by magnetic reconnection or bouncing in the foot point. These are intriguing 
features which may characterize the foot points, whereas their energy is not dominant. Such 
features are not predicted by numerical simulations (Takahashi et al. 2009) and more realistic 
simulations are desirable to pursue these details. 

6. Conclusions 

We have carried out 12 CO(J=2-l) and 12 CO(J=3-2) observations at spatial resolutions 
of 1.0-3.8 pc toward the entirety of loops 1 and 2 and part of loop 3 in the Galactic center 
with NANTEN2 and ASTE. These new results revealed detailed distributions of the molecular 
gas and the line intensity ratio of the two transitions, #3-2/2-1- Main points of the study are 
summarized below; 

1. The 12 CO J =3-2 and 2-1 transitions show similar distributions to 12 CO J =1-0 in loops 
1 and 2. In the three loops, the line intensity ratio #3-2/2-1 is in a range from 0.1 to 2.5 
with a peak at ~0.7. The ratio in the disk molecular gas is in a range from 0.1 to 1.2 with 
a peak at 0.4, significantly smaller than that in the Galactic center. This suggests that 
the loops are more highly excited than the disk molecular gas. 

2. An LVG analysis of three transitions, 12 CO J =3-2 and 2-1 and 13 CO J =2-1, toward six 
positions in loops 1 and 2 shows that the density and temperature are in a range 10 2,2 - 
10 4 ' 7 cm -3 and 15-100 K or higher, respectively. Considering the frequency distribution 
of #3-2/2-1, we suggest that the loops are characterized by higher temperature and/or 
density significantly different from those in the Galactic disk molecular gas. 

3. We have found that three regions extended by 50-100 pc in the loops tend to have higher 
excitation conditions as characterized by #3-2/2-1 greater than 1.2. These include a region 
inside loop 1 (Region B), a region in the foot points of loops 1 and 2 (Region C) and a region 
above the foot point of loop 2 (Region D). In addition, we find two hot compact spots with 
high # 3 _2/2-i above Region C. The highest ratio of 2.5 is found in the most developed foot 
points in Region C. The high ratio in the foot points is interpreted as evidence for strongly 
shocked conditions which is consistent with their large linewidths of 50-100 km s _1 . The 
other two regions outside the foot points suggest that the molecular gas is heated up by 
some additional heating mechanism other than the gravitational energy release possibly 
including magnetic reconnection, since the regions have no associated young massive stars. 

4. A detailed analysis of four of the foot points have shown a U shape, an L shape or a 
mirrored-L shape in the latitude-velocity distribution. These shapes share similarity with 
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the U shape reported in the foot points between loops 1 and 2 by TlOb. The U and 
L/mirrored-L shaped features have large velocity widths of 20-100 km s _1 and have broad 
features in their bottom. It is shown that a simple kinematical model which incorporates 
global rotation and expansion of the loops in addition to falling motion along the loops is 
able to explain the characteristic shapes. 
5. We found that the foot points have localized broad features that appear to lie between 
the two velocity features constituting the U shapes. Because these broad features have 
no associated young stars, they require some non stellar energy source, where the typical 
kinetic energy involved in one of such broad features is 10 49 erg. We argue that the 
broad compact features may be due to magnetic reconnection in the U shaped gas where 
magnetic field lines in the opposite directions become reconnected. We present also found 
protrusions in the bottom of foot points possibly formed by magnetic reconnection or 
bouncing in the foot points. 
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Appendix 1. Velocity channel distributions of 12 CO(J =2-1), 12 CO(J =3-2) and 
-R3-2/2-1 

Figures A and B show velocity channel distributions of 12 C0( J =2-1) and 12 C0( J =3- 
2) in I and b, respectively 2 . Loop 1 is for a velocity range from —160 km s _1 to —50 km s" 1 , 
loop 2 from —120 km s" 1 to —30 km s _1 and loop 3 from km s" 1 to 120 km s _1 . The thick 
contours show the lowest 3a level of the 12 CO( J =1-0) intensity, 5.0 K km s _1 . It is not shown 
in the present paper but we note that 12 CO(J =3-2) emission is detected in loop 3 top in 
(l,b) ~(357.°5-356.°8, 0.°3-l.°7) and v =60-120 km s _1 . Loops 1 and 2 tend to become separated 
from the center with the increase of velocity while loop 3 shows an opposite trend. Fujishita 
et al. (2009), and TlOa give a simple kinematical model to explain this velocity distribution by 
assuming rotational motion with expansion. 

The local features outside the Galactic center have lower velocities of —40-20 km s" 1 
and they show different distributions in 12 CO(J=1-0), 12 CO(J=2-l) and 12 CO(J=3-2). The 
b distribution of the 12 CO(J =1-0) and 12 CO(J =2-1) is extended to b =1.0°, whereas the 
12 CO(J =3-2) is limited to b less than 0.°25 in 0-20 km s _1 . This difference suggests that 
the 12 CO (J =3-2) selectively traces denser/warmer gas near the plane than the lower J lines 
and we may separate the loop features from the foreground local features at higher angular 
resolutions. 

Figure C shows velocity channel distributions of the line intensity ratio #3-2/2-1 for loops 

1 and 2. The regions where #3-2/2-1 is greater than 1.5 are four regions; 1) toward the loop 1 
top in v ~ —80 — 70 km s~ x and (l,b) ~(356.°5, l.°2), 2) toward the loop 2 top in v ~ —60 — 50 
km s" 1 and (l,b) ~(355.°3, l.°5), 3) toward the loop 1 west top in (l,b) ~(356.°2, 0.°9), and 4) in 
v ~ 0-10 km s" 1 around (l,b) ^(355.°2, l.°5). The fourth region has very small linewidths of a 
few km s _1 in 12 CO( J =3-2) in the velocity range of local features and we shall not deal with 
this in the present work. 

Appendix 2. Modeling the velocity distribution of the foot points 

We make a simple kinematic modeling of the three loops. 

Figure D presents the model for loops 1 and 2 where we adopt a radius of the loops 670 pc 
with an expansion velocity of 150 km s _1 and a rotational velocity of 50 km s" 1 (Section 6.3 of 
TlOa). The position angle of 55° is adopted for the position angle between loops (Figures D(a) 
and D(b)). The falling velocity and the loop shape are taken from the numerical simulation by 
Takahashi et al. (2009) at an elapsed time of 5.25 x 10 7 yr so as to fit the loop height (Figure 
D(b)). A l-b diagram and a b-v diagram are shown in Figures D(d) and D(e). Figure E is for 
loop 3 whose radius is 1000 pc (Section 5.4 of Fujishita et al. (2009)) with a position angle of 

2 The fits files of loops 1 and 2 in 12 CO(J=2-l) and 12 CO(J=3-2) can be taken from 
http : / / www. a. phys . nagoya-u. ac . j p/ ~ torii/loop/cg 
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the center of the loop, 145° (Figures E(a) and E(b)). The expansion velocity is 130 km s _1 and 
the rotational velocity is 80 km s -1 . The falling velocity and the shape is tentatively taken from 
a model at 7.84 x 10 7 yr from Takahashi et al. (2009) to fit the loop height (Figure E(b)). An 
l-b diagram and a b-v diagram are shown in Figures Ed and E(e). The dotted lines in Figures 
D and E(e) indicate that without falling motion. 

These model calculations show that the foot points show U or L/mirrored-L shaped 
distributions in b-v diagrams, independently of the falling motion. This indicates that the U 
and L/mirrored-L shapes basically reflect expansion and rotation. Figures D and E(e) show 
boxes of Regions A, C, E and G. The trend of U and L/mirrored-L shapes is qualitatively 
reproduced as is consistent with the observations. We, however, note that the quantitative 
agreement is not currently good. This is perhaps because the model by Takahashi et al. (2009) 
has no fine tuning to fit the observations and, rather, intended to reproduce physical processes 
on the loops. 
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Table 1: Positions of hot regions in -90 - -40 km s 



ID b( )* maxi? 3 _2/2-i reference 1 ' 



p 


356.40 


1.24 


1.90 


Region B 


Q 


356.28 


0.92 


1.28 


Region C 


R 


356.04 


1.54 


1.55 


Region C 


S 


355.89 


1.90 


1.15 


Region D 


T 


355.37 


1.76 


1.27 


Region D 


W 


355.23 


1.38 


1.53 


Region D 



* These positions indicate maximum -K3-2/2-1 positions. 
t These regions are shown in Figure 4. 



Table 2: Positions of 13 CO( J = 2-1) observations 



ID 


I 


b 




12 CO(J 


= 2-1) 




(°) 


(°) 


Tpeak (K) 


V P eafe(km s _1 ) 


AV(km s- 1 ) 


Hot regions 


a 


357.58 


0.42 


1.7 


-139.5 


14.0 




1) 


356.55 


1.33 


3.2 


-67.3 


39.7 


P 


cl 


356.25 


0.84 


7.0 


-90.0 


58.0 


Q 


c2 


356.06 


1.52 


0.8 


-94.1 


48.3 


R 


cl 


355.38 


1.57 


1.9 


-62.0 


26.5 


U 


c 


355.53 


0.80 


3.9 


-75.7 


53.5 
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Table 3: i?3_2/2— l and -R13/12 in peaks a-e 



Peak 


Vlsr range 


Line ratio 


(km s- 1 ) 


-R3-2/2-1 


-R13/12 




-155 — 145 


0.52 


0.05 




-145 — 135 


1.05 


0.08 


a 


-135 — 125 
-65—55 


0.85 
0.35 


0.06 
0.05 




-45—35 


0.53 


0.04 




-35 — 25 


0.38 


0.02 


1) 


-75—65 
-65 — 55 


0.72 
1.03 


0.08 
0.09 




-90 — 80 


0.97 


0.11 




-80 — 70 


1.20 


0.10 


cl 


-70 — 60 
-60 — 50 


0.92 
0.60 


0.10 
0.06 




-40 — 30 


0.33 


0.06 




-30 — 20 


0.37 


0.09 




-100 — 90 


1.10 


0.37 


c2 


-90 — 80 


0.80 


0.19 




-80 — 70 


0.56 


0.14 


d 


-70 — 60 
-60 — 50 


1.18 
1.50 


0.05 
0.07 




-100 — 90 


0.92 


0.03 




-90 — 80 


0.95 


0.04 




-80 — 70 


0.80 


0.03 




-70 — 60 


0.63 


0.04 


c 


-60 — 50 


0.53 


0.03 




-50 — 40 


0.38 


0.06 




-40 — 30 


0.33 


0.02 




0-10 


0.35 


0.05 




10-20 


0.62 


0.05 



* The averaged values in a 10 km s 1 winows. 
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Table 4: LVG results for X(CO) / (dv / dr) = 1.1 xlCT 5 pc (km s^ 1 )" 1 



Peak 


Vlsr 


n(H 2 ) (cm 


- 3 ) 


T k 


(K) 


2 


(km s _1 ) 


X 2 < 6.25 


min. x 2 


X 2 < 6.25 


min. x 2 


min 




-155 — 145 


10 2.7_ 10 3.i 


10 2.9 


19-95 


42 


0.4 x 10~ 2 




-145 — 135 


1Q 3.6 < 


1Q 4.4 


101 < 


195 


1.7 x 10~ 2 


M 

a 


-135 — 125 


10 3.0_ 10 3.7 


10 3.2 


47-249 


114 


0.2 x 10~ 2 


-65 — 55 


10 2.8_ 10 3.0 


1Q 2.9 


12-34 


18 


0.6 x IO" 2 




-45 — 35 


10 2.7_ 10 3.0 


10 2.8 


34-127 


60 


0.07 x 10~ 2 




-35 — 25 


< io 2 - 3 


— 


123 < 


— 


2.1 


h 


-75 — 65 


10 3.1_ 1Q 3.5 


10 3.2 


25-86 


44 


0.4 x 10~ 2 


l J 


-65 — 55 


1Q 3.5 < 


io 4 - 7 


70 < 


214 


1 x 10~ 5 




-100 — 90 


— 


— 


— 


— 


31.2 




-90 — 80 





— 


— 


— 


5.8 




-80—70 


10 2.8„ 10 3.1 


10 2.9 


167 < 


— 


3.4 x 10~ 2 




-70—60 


10 2.6_ 10 2.9 


io 2 - 7 


87 < 


194 


0.3 x 10" 2 


cl 


-60 — 50 


1 2.4_ 10 2.8 


10 2.6 


53 < 


128 


0.01 x 10~ 2 




-50 — 40 


10 2.0_ 10 2.7 


10 2.2 


37 < 


188 


1 x IO" 5 




-40 — 30 


< io 2 - 6 


— 


31 < 


245 


1.3 




0-10 


10 2.8_ 10 3.0 


1Q 2.9 


12-26 


16 


1.7 x 10~ 2 




10-20 


10 2.6_ 1Q 3.6 


10 30 


32-255 


48 


0.05 x 10~ 2 




-100 — 90 





— 




— 


60.0 


c2 


-90 — 80 





— 


— 


— 


7.0 




-80—70 


io 3 - 4 < 


10 3.8 


21-190 


66 


0.2 x 10~ 2 


u 


-70 — 60 
-60 — 50 


10 3.4_ 10 4.5 




— 
— 


245 < 
— 


— 
— 


0.4 
14.5 




-90—80 


1Q 3.5 < 


1Q 3.9 


45-280 


94 


0.08 x 10~ 2 




-80—70 










5.8 




-70—60 


io 3 - 4 < 


1Q 3.7 


45-286 


82 


0.04 x 10~ 2 


e 


-60 — 50 


10 3 - io 3 - 2 


10 31 


22-71 


36 


0.03 x 10~ 2 




-40 — 30 


10 2.8 _ 1Q 3.1 


1Q 3.0 


< 24 


14 


1.8 x IO" 2 




-30—20 


1Q 3.1 _ 1Q 3.3 


10 3.2 


< 15 




0.3 x IO" 2 
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Table 5: Physical properties of the compact broad components in panel b of Figures 9, 10, 12 and 13 



Region 


no. 


^CO(J=3-2) (10 3 


K km S 1 ) 


°3-2/2-l 


^3-2/1-0 


Mass (10 3 M S ) t 


A 


1 

2 


0.5 
6.6 




0.23 
0.40 







— 
— 




3 


3.5 




0.77 





— 


c 


4 
5 


0.4 
0.6 




* 

0.93 


0.58 
0.46 


2.0 
1.5 




6 


2.2 




0.90 


— 


— 




7 


5.8 




— 


— 


— 


E 


8 


2.1 




0.80 


— 


— 




9 


4.1 




0.73 


— 


— 




10 


4.8 






— 


— 




1 1 

11 


4.0 




0.50 






G 


12 


0.7 




0.55 








13 


0.7 




0.52 








14 


1.4 




0.55 






Notes. '-' 
* 12 CO(J 


means 
=2-1) 


no value because ot outside ot the obscr 
is not detected at more than 3<r level. 


ving region in the 


line. 





t Using Mopra 12 CO(J=1-0) data (TlOb). 
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Color, Thin contours 12 CO(J=2-1) 

n '_ Thick contours 12 CO(J=1-0) J 

U D i i i , i , , i , i , i i , 

358.0 357.0 356.0 



Galactic Longitude [degree] 



Fig. 1: 12 CO(J=2-l) integrated intensity distribution of loops 1 and 2. (a)Integrated intensity 
distribution of loop 1. The integration range in velocity is from —180 km s~ x to —90 km s" 1 . 
(b)Integrated intensity distribution of loop 2. The integration range in velocity is from —90 
to —40 km s _1 . In both figures, color and thin contours show the 12 CO(J =1-0) emission and 
thick contours in black and gray show the 12 CO( J =1-0) emission. Here, thin contours in black 
show the loops and in gray show the other components. Thin contours are plotted every 18 K 
km s~ l from 4.8 K km s _1 , and thick contours are plotted at 7 K km s _1 and 100 K km s" 1 . 
Dashed lines indicate the observing region of 12 CO(J =2-1). 
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Fig. 2: 12 CO(J =3-2) integrated intensity distribution of loops 1 and 2. (a)Integrated inten- 
sity distribution of loop 1. The integration range in velocity is from —180 to —90 km s _1 . 
(b)Integrated intensity map of loop 2. The integration range in velocity is from —90 to —40 
km s" 1 . In both figures, color and thin contours show the 12 CO(J =1-0) emission and thick 
contours in black and gray show the 12 CO( J =1-0) emission. Here, thin contours in black show 
the loops and in gray show the other components. Thin contours are plotted every 30 K km 
s _1 from 8 K km s _1 , and thick contours are plotted at 7 K km s _1 and 100 K km s _1 . Dashed 
lines indicate the observing regions of 12 CO(J =3-2). 
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Fig. 3: CO integrated Intensity distributions of loop 3. The integration range in velocity is 
from 30 to 160 km s" 1 . Color and thin contours in (a) show the 12 CO(J =2-1) emission and 
those in (b) the 12 CO(J =3-2) emission. The levels of the contours are the same as that in 
Figures 1 and 2. In both figures, thick contours in black and gray show the 12 CO(J =1-0) 
emission. Here, thin contours in black show the loops and in gray show the other components. 
Dashed lines in (a) and (b) show the observing region of 12 CO(J =2-1) and 12 CO(J =3-2), 
respectively. 
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Color: 12 CO(J=3-2), Contours: 12 CO(J=1-0) C 
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Fig. 4: 7 regions from A to G defined in the text are shown by the boxes on the CO integrated 
intensity distributions of loops 1, 2 and 3. Color and contours show the integrated intensity 
distributions of 12 CO(J=3-2) and 12 CO(J=1-0), respectively 
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Color R3-2/2-1, Contour 12 CO(J=1-0) 




358.0 357.0 356.0 

Galactic Longitude [degree] 



Fig. 5: Distributions of the integrated intensity ratio from 12 CO(J =2-1) to 12 CO(J =3-2), 
i?3_ 2 /2-i- (a)i?3-2/2-i distributions of loop 1. The integration range in velocity is from —180 to 
—90 km s -1 . (b)i?3_2/2-i distributions of loop 2. The integration range in velocity from —90 to 
—40 km s _1 . Contours in both figures show the 12 CO(J=1-0) integrated intensity distribution 
and are plotted at 7 K km s" 1 . Dashed lines indicate the 12 CO( J =2-1) observing region. 
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Fig. 6: i?3_2/2-i distributions of loop 3. The integration range in velocity is from —180 to —90 
km s _1 . Contours show the 12 CO(J =1-0) integrated intensity distribution and are plotted at 
7 K km s _1 . Dashed lines indicate the 12 CO(J=2-l) observing region. 



30 




Fig. 7: Intensity weighted frequency distributions of R3-2/2-1 in loops 1 and 2. The data are 
smoothed to a 90" spatial resolution with a gaussian function and smoothed to a 2 km s" 1 
velocity resolution. The intensity threshold of 12 CO(J =2-1) and 12 CO(J =3-2) is 3cr. The 
white lines show the contributions of the loops ( — 180 km _1 to —40 km ), and the red, orange 
and green lines show the contributions of the local components with a interval of 20 km _1 from 
—40 km -1 . The high -R3— 2/2— 1 region between km _1 and 15 km _1 was subtracted from the 
results (see Figure C in Appendix 1 of Supporting Online Material). 
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(a) loops 1 & 2 



loops 1&2 all 
(-180 --40 km s' 1 ); 
Region A 
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1 .0 2.0 

12 CO(J=3-2)/ 12 CO(J=2-1) 
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Region F 
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1.0 2.0 
12 CO(J=3-2)/ 12 CO(^2-1) 



3.0 



Fig. 8: Intensity weighted frequency distributions of i?3_2/2-i- The data are smoothed to a 
90" spatial resolution with a gaussian function and smoothed to a 2 km _1 velocity resolution. 
The intensity threshold of 12 CO( J =2-1) and 12 CO( J =3-2) is 3a. Colored lines in the figures 
show the contributions of Regions A-G shown in Figure 7. (a) i? 3 _ 2 /2-i frequency distributions 

of loops 1 and 2. The filled black region show the whole of loop 1 ( — 140 40 km _1 ). Red 

shows Region A (the east footpoint of loop 1), pink shows Region C (the right footpoint of loop 
1) at b < 1.5°, yellow shows Region C at b > 1.5°, green shows Region D (the top of loop 2), and 
blue shows Region E (the right footpoint of loop 2). (b) #3-2/2-1 frequency distributions of 
loop 3. The filled black region shows the whole of loop 3 (30 - 160 km ~ 1 ). Red shows Region 
F (the loop 3 west footpoint) and orange shows Region G (the top of loop 3). 
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Fig. 9: (a-top) Velocity channel maps of 12 CO(J=3-2) integrated intensity every 10 km s _1 
toward the east footpoint of loop 1. Contours are illustrated every 10 K km s~ x from 4 K km 
s" 1 . (a-bottom) Velocity channel distributions of Rz-2/2-1 averaged every 10 km s _1 . Contours 
show the 12 CO(J=3-2) emission smoothed to a 90" spatial resolution with a gaussian function 
and are plotted every 10 K km s _1 from 4 K km s _1 . (b-left) Integrated intensity distribution 
of 12 CO(J=3-2) in Region A. An X-Y coordinate system is defined in the text. Inclination of 
X axis is ~45°. (b-right) X axis channel map of the 12 CO(J=3-2) in region A . Contours are 
plotted every 0.05 K deg. from 0.02 K deg. Hereafter, dashed-lines indicated the 12 CO(J=3-2) 
observing regions. 
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Fig. 10: (a-top)Velocity channel distributions of 12 CO(J=3-2) integrated intensity every 10 
km s" 1 at the west footpoint of loop 1. Contours are illustrated every 20 K km s _1 from 6 K 
km s~ x . (a-bottom) Velocity channel distributionss of R3-2/2-1 integrated integrated intensity 
every 10 km s _1 . Contours show the 12 CO(J=3-2) emission smoothed with a gaussian function 
to a 90" spatial resolution and are plotted every 20 K km s _1 from 6 K km s _1 . 
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Fig. 10: (b)Longitude velocity distributions of the 12 CO(J=3-2) with interval of 100" in I in 
regions C. Contours are plotted every 0.05 K deg. from 0.02 K deg.. 
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Color, Thin contour : 12 CO(J=3-2) 50pc 
Thick contour : 12 CO( J=1 -0) j , , H 




Color : R3-2/2-1, Contour : 1g CO(J=3-2) (9CT smoothed) 




Galactic Longitude [degree] 

Fig. 11: (top) Integrated Intensity distributions of loop 2 in 12 CO( J=3-2), color and thin 
contours, and 12 CO(J=1-0), thick contours, in region D. The integration range in velocity 
is from —90 km s _1 to —60 km s" 1 . Dashed lines indicate the observed regions, (bottom) 
Distributions of the #3-2/2-1 ■ Contours show the 12 CO(J=3-2) emission smoothed with a 
gaussian function to a 90" spatial resolution and are plotted every 10 K km s" 1 4 K km s^ 1 . 
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Fig. 12: (a-top) Velocity channel distributions of 12 CO(J=3-2) integrated intensity every 10 
km s _1 at the west footpoint of loop 2. Contours are illustrated every 10 K km s _1 from 4 K 
km s _1 . (a-bottom) Velocity channel distributions of the R3-2/2-1 integrated every 10 km s _1 . 
Contours show the 12 CO(J=3-2) emission smoothed with a gaussian function to a 90" spatial 
resolution and are plotted every 10 K km s~ x 4 K km s _1 . (b)Longitude velocity distributions 
of the 12 CO(J=3-2) with interval of 100" in I in regions E. Contours are plotted every 0.05 K 
deg. from 0.02 K deg.. 
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Fig. 13: (top) Velocity channel distributions of 12 CO(J=3-2) integrated intensity every 10 km 
s" 1 at the west footpoint of loop 3. Contours are illustrated every 10 K km s" 1 from 4 K km s" 1 . 
(bottom) Velocity channel distributions of the #3-2/2-1 intensity ratio integrated every 10 km 
s _1 . Contours show the 12 CO(J=3-2) emission that was smoothed with a Gaussian function 
to a 90h spatial resolution and are plotted every 10 K km s" 1 from 4 K km s" 1 . 
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g. 13: (b)Longitude velocity distributions of the 12 CO(J=3-2) with interval of 100" in I 
gions G. Contours are plotted every 0.05 K deg. from 0.02 K deg.. 
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. Bold line : X/(dv/dv) = 1 .1 xl0" 5 pc (km s" 1 )" 1 
. Dashed line : X/(dv/dr) = 1 .1 xl0" 6 pc (km s" 1 )" 1 




Fig. 14: i? 3 _2/2-i as a function of kinetic temperature estimated by the LVG calculations. 
We take two X(CO) / \dv / dr) , 1.1 x l(T 5 (Bold line) and 1.1 x 10" 6 (Dashed line). Colored lines 
show densities of 10 2 '°(green), 10 2,5 (blue), 10 3 ' (purple), 10 3,5 (red) and 10 4,0 cm~ 3 (black). 
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Fig. 15: CO spectra at 6 positions shown in Figure 5. The intensities of the 13 CO( J =2-1) 
emission were multiplied by 5. 
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Fig. 16: Chi square, x 2 ; distributions estimated from the LVG calculations on the temperature- 
density plane for the different velocity of the spectra shown in Figure 15. Colored lines are 
plotted at x 2 of 3.84, which corresponds to a 95% confidence level, and the crosses show the 
points that minimum \ 2 is found. 
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Fig. 17: LVG results for X(CO) / \dv / dr) of 1.1 x 10 -5 for six spectra shown in Figure 15. 
Horizontal axis of all figures is velocity. Top row shows the 12 CO(J=3-2) spectra at the peaks 
a-f shown in Figure 5. The second and third rows show number density, n(B.2), and kinetic 
temperature, Tjt, respectively. Open circles show the lowest point of x 2 i and the error range is 
defined as 5% confidence level of x 2 distribution with 1 degree of freedom, which corresponds 
to x 2 = 3.84. 
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Fig. 18: The sets of the longitude- latitude distributions (left panels) and velocity-latitude 
distributions (right panels) in 12 CO( J =3-2) at the four footpoints of the loops. In the Region 
A, the right panel is the velocity distributions for the X-axis which was defined in the left panel 
of the figure. 
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Fig. 19: (left) A schematic image of the positions of loops 1, 2 and 3 in a face on view, (right) 
Latitude - velocity diagram of the four foot points, top-panels: Outer boundaries of the foot 
points are shown. Origins of the perpendicular axis of the regions A, C, E and G are b =0.°68, 
0.°8, 0.°8 and 0.°56, respectively, bottom-panels: Schematic images of U shape. 
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Fig. A: Velocity-channel distributions of the 12 CO(J=2-l) emission integrated every 10 km s _1 
from —180 to 60 km s _1 . Color and thin contours (black and white) show the 12 CO(J =2-1) 
emission. The contours are plotted every 6 K km s -1 (black) and 12 K km s" 1 (white) from 4.7 
K km s _1 . Thick contours show the 12 CO (J =1-0) emission and are plotted at 7 K km s _1 . 
The numbers shown in the left-top of the each panel show the integration range in km s _1 . 
Dashed-lines indicate the observed region. <r 
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Fig. A: Continued, but the levels of thin contours are different from the previous figure. The 
contours are plotted every 12 K km s _1 (black) and 24 K km s~ 1 (white) from 4.7 K km s" 1 . 
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Fig. A: Continued. 
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Fig. A: Velocity-channel distributions in 12 CO(J =2-1). Thin black and white contours are 
12 CO(J =2-1), plotted every 24 K km s _1 (black) and 48 K km s _1 (white) from 4.7 K km s" 1 . 
Dashed-lines indicated the observed region. 
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Thick contour : 12 CO(J=1-0) ASTE NANTEN 100pc 25 50 
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Fig. B: Velocity-channel distributions of the 12 CO(J=3-2) emission integrated every 10 km s _1 
from —180 to 120 km s" 1 . Color and thin contours show the 12 CO(J =3-2) emission, and the 
contours are plotted every 10 K km s _1 from 4 K km s _1 . Thick contours show the 12 CO(J=l- 
0) emission and are plotted 7 K km s _1 . The numbers shown in the left-top of the each panel 
show the integration range in km s -1 . Dashed-jlines indicate the observed regions. 
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Fig. B: Continued, but the levels of the thin contours are different from the previous figure, 
that are plotted every 20 K km s _1 from 4 K km s" 1 . 
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Fig. B: Continued. 
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Fig. B: Color and thin contours are 12 CO( J =3-2), plotted every 10 K km s 1 from 4 K km s 1 
and thick contours are 12 CO (J =1-0), plotted 7 K km s _1 . Dashed-lines indicate the observed 



regions. 
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Fig. B: Continued. 
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Fig. C: Velocity channel distributions of R3-2/2-1 integrated every 10 km s _1 . Thin contours 
show the 12 CO(J =2-1) emission and are plotted at 4.7 K km s" 1 . Thick contours show the 
12 CO(J =1-0) emission and are plotted at 7 K km s" 1 . Dashed-lines indicate the analysed 
regions. 
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Fig. C: Continued. 
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Fig. C: Continued. 
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Fig. C: Continued. 
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Fig. D: (a)Schematic image of the distribution of loops 1 and 2 in a face-on view. (b)Model of 
two loops on the cylinder coordinate. (c)Down flow speed parallel to the galactic plane on the 
loops calculated by Takahashi et al. (2009). (d)Model calculation of the Galactic Longitude- 
Latitude diagram of two loops estimated with the model, (e) Model calculation of the latitude- 
velocity diagram of the loops. Dashed-lines indicate the result without down flow speed. Boxes 
indicate the footpoints of loops 1 and 2 (Regions A, C and E). 
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Fig. E: (a) Schematic image of the distribution of loop 3 in a face on view, (b) Model of a loop 
on the cylinder coordinate. (c)Down flow speed parallel to the galactic plane following the 
loop calculated by Takahashi et al. (2009). (d)Model calculations of the Galactic longitude- 
Latitude diagram of the loop. (e)Model calculation of the latitude-velocity diagram of two 
loops. Dashed- lines indicate the result without down flow speed. Boxes indicate the footpoints 
of loop 3 (Region G). 
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